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We consider the interpretation of the MiniBooNE low-energy anomaly and the Gallium radioac- 
tive source experiments anomaly in terms of short-baseline electron neutrino disappearance in the 
framework of 3-1-1 four-neutrino mixing schemes. The separate fits of MiniBooNE and Gallium data 
are highly compatible, with close best-fit values of the effective oscillation parameters Am^ and 
sin^2i?. The combined fit gives Am^ > 0.1 eV^ and 0.11 < sin^ 2i9 < 0.48 at 2a. We consider 
also the data of the Bugey and Chooz reactor antineutrino oscillation experiments and the limits 
on the effective electron antineutrino mass in /3-decay obtained in the Mainz and Troitsk Tritium 
experiments. The fit of the data of these experiments limits the value of sin^ 2i9 below 0.10 at 
2a. Considering the tension between the neutrino MiniBooNE and Gallium data and the antineu- 
trino reactor and Tritium data as a statistical fluctuation, we perform a combined fit which gives 
Am^ ~ 2eV and 0.01 < sin^ 2i? < 0.13 at 2a. Assuming a hierarchy of masses mi, m2, m,3 ^ m4, the 
predicted contributions of m^ to the effective neutrino masses in /3-decay and neutrinoless double- 
/3-decay are, respectively, between about 0.06 and 0.49 and between about 0.003 and 0.07 eV at 2a. 
We also consider the possibility of reconciling the tension between the neutrino MiniBooNE and 
Gallium data and the antineutrino reactor and Tritium data with different mixings in the neutrino 
and antineutrino sectors. We find a 2.6a indication of a mixing angle asymmetry. 

PACS numbers: 14.60. Pq, 14.60.Lm, 14.60. St 
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I. INTRODUCTION 

Neutrino oscillations have been observed in solar, at- 
mospheric and long-baseline reactor and accelerator ex- 
periments. The data of these experiments are well fitted 
in the framework of three-neutrino mixing, in which the 
three flavor neutrinos j/g, v^, Vt are unitary linear com- 
binations of three massive neutrinos vi, 1/2, 1^3 with the 
solar (SOL) and atmospheric (ATM) squared-mass dif- 
ferences 



Am^i = Ato^ol - 8 X 10"^ eV^ 



lAr 



'32 



= Ar 



'ATM 



2 X 10"^ cV' 



(1) 
(2) 



m^ and rrij is the mass of the neu- 



where AmL 

trino Vj (see Refs. [lH8|)- 

Besides these well-established observations of neutrino 
oscillations, there are at least three anomalies which 
could be signals of short-baseline neutrino oscillations 
generated by a larger squared-mass difference: the LSND 
Dfj^ — >■ De signal [91, t he Gallium radioactive source exper- 
iments anomaly [Ifl, [ill , and the MiniBooNE low-energy 
anomaly [l3]. In this paper we consider the MiniBooNE 
and Gallium anomalies, which can be explained by short- 
baseline electron neutrino disappearance [13h15i] in the 
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effective framework of four-neutrino mixing, as explained 
in Sections M and Uni On the other hand, the LSND 
anomaly is disfavored by the results of the MiniBooNE 
i/fj^ — ^ fe ex per iment [121 [l6j and may require another 
explanatio n ll7H23j . 

In Refs. [ij,[la| we propo sed to explain the MiniBooNE 
low-energy anomaly |12l . ll6l | through the disappearance of 
electron neutrinos due to very-short-baseline oscillations 
into sterile neutrinos generated by a squared-mass differ- 
ence Am^ larger than about 20 eV^. In that case, the 
analysis of the MiniBooNE data is simplified by the fact 
that the effective survival probability P^^^^^ is practi- 
cally constant in the MiniBooNE energy range from 200 
to 3000 MeV. In this paper we extend the analysis of 
MiniBooNE data to lower values of Am^, considering 
the resulting energy dependence of the effective short- 
baseline (SBL) electron neutrino and antineutrino sur- 
vival probability 



IlTaL^E) 



SVC? 2-0 sin^ 



Am^L 
AE 



(3) 



where L is the neutrino path length and E is the neu- 
trino energy (OPT invariance implies that the survival 
probabilities of neutrinos and antineutrinos are equal; see 
Ref. i). 

The two-neutrino-like effective short-baseline survival 
probability in Eq. ^ is obtained in four-neutrino 
schemes (see Refs. [H, 0, B 0]); which are the simplest 
extension of three- neutrino mixing schemes which can ac- 



commodate the two small solar and atm.ospheric squared- 
mass differences in Eqs. ([T|) and ^, and one larger 
squared-mass difference for short-baseline neutrino oscil- 
lations, 



lAm^il = Am^ >0.1eV^ 



(4) 



The existence of a fourth massive neutrino corresponds, 
in the flavor basis, to the existence of a sterile neutrino 

In this paper we consider 3+1 four-neutrino schemes, 
since 24-2 four-neutrino schemes are disfavored by the 
combined constraints on active-sterile transitions in solar 
and atmospheric neutrino experiments [J] . For simplicity, 
we consider only 3+1 four-neutrino schemes with 



mi, 7712, ms ^ rn4 . 



(5) 



which give the Arrili in Eq. ((4]) and appear to be 
more natural than the other possible 3+1 four-neutrino 
schemes in which either three neutrinos or all four neu- 
trinos are almost degenerate at a mass scale larger than 



VAm2 (see Refs. [111 SB) 

In 3+1 four- neutrino schemes the effective mixing an- 
gle in the effective short-baseline electron neutrino sur- 
vival probability in Eq. ^ is given by (see Refs. [l|,|j, @, 

H) 



sin2 2l? = 4|C/e4p(l-|C/e4| 



(6) 



In this paper we assume that the value of |L''^4p is so 
small that the effective short-baseline muon neutrino sur- 
vival probability is practically equal to unity and short- 
baseline I'f^ ^ i^e are negligiblqj. This assumption is jus- 
tified by the lack of any indication of f^ — > i>e transitions 
in the MiniBooNE experiment [l^, UM ^^'^ the limits on 
short-baseline muon neutrino disappearance found in the 
CDHSW H, CCFR [H and MiniBooNE f2| experi- 
ments. We do not consider the MiniBooNE antineutrino 
data |27, ], which have at present statistical uncertainties 
which are too large to constraint new physics [15| . 

The plan of the paper is as follows. In Section |lT] we 
discuss the analysis of MiniBooNE data. In Section IIIII 
we present an update of the analysis of Gallium data 
published in Ref. IJ] and the combined analysis of Mini- 
BooNE and Gallium data. In Section ITVl we discuss the 
implications of the measurements of the effective electron 
neutrino mass in Tritium /?-decay experiments and their 
combination with reactor neutrino oscillation data. In 



^ In 3+1 four-neutrino schemes the effective short-baseUne muon 
neutrino survival probability has the form in Eq. ^ with 
sin^ 2i9 replaced by sin^ 2i?^^ = 4\U ^4]'^ {l - \U ^4]^) . The ef- 

transition probability is given 



(-) 



fective short-baseline fj 

sin^ 2»9, 



MM 

(- 



4|f/e4p|C/M4p (see Refs. [E1IS0) 



2/ Am/L' 



with sin^ 2i?e 



Section |V] we present the results of the combined anal- 
ysis of MiniBooNE, Gallium, reactor and Tritium data 
and in Section I VI I we present the corresponding predic- 
tions for the effective masses measured in /3-decay and 
neutrinoless double- /3-decay experiments. In Section IVlII 
we calculate the mixing angle asymmetry between the 
neutrino and antineutrino sectors which could explain 
the tension between the neutrino and antineutrino data 
under our short-baseline i^g-disappearance hypothesis. In 
Section rVlIII we draw the conclusions. 



II. MINIBOONE 



The MiniBooNE experiment was made with the pur- 



pose of checking the indication of Vf^ 



oscillations 



generated by a Am^ > 0.1 eV^ found in the LSND ex- 
periment . The MiniBooNE collaboration did not find 
any indication of such oscillations in the Vfj, — )■ Vg chan- 
nel [11,13. On the other hand, the MiniBooNE coUab- 
oration found an anomalous excess of low-energy i^e-like 



events in the data on the search for 1/^, 



oscillations 



[H, ill], as shown in Fig. [T^. 

As in Refs. jTS, llS], we consider an explanation of the 
low-energy MiniBooNE anomaly based on the possible 
short-baseline disappearance of electron neutrinos, tak- 
ing into account a possible overall normalization factor fi, 
of the calculated ^'e-induced and misidentified ^'^-induced 
events which contribute to the observed number of j/g- 
like events. The normalization factor f^, could be due 
mainly to the uncertainty of the calculated neutrino fiux 
(see Ref. [2a| ) . Since the misidentified t'^-induced and j/g- 
induced events dominate, respectively, at low and high 
energies (see Fig. [T^), the low-energy excess can be fit- 
ted with /^ > 1 and the high-energy data can be fitted 
compensating f^, > 1 with the disappearance of I'e's. 

In Refs. [l^iia] we considered only very-short-baseline 
i^e disappearance due to a Arn^ > 20 eV^, which gener- 
ates a survival probability P,^^^^^ which is constant in 
the MiniBooNE energy range, from 200 to 3000 MeV. 
In this paper we extend the analysis to lower values of 
Am^, considering the resulting energy dependence of the 
survival probability. In this case, the theoretical number 
of j/g-like events in the jth energy bin is given by 



AT-thc jyrl^e,iiiC . ^I',, ,thG 



where 



j^u,,thc ^ f pU) jv'^^' 



cal 

i 



is the number of z^„-induced events and 






(7) 



(8) 



(9) 



is the number of misidentified j/^-induced events 



Nf^f" and iV." 



cal 



Here 
are, respectively, the number of v^- 



induced and misidentified i/^-induced events calculated 
by the MiniBooNE collaboration for the jth energy bin 
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FIG. 1. Expected number of Ue events compared with MiniBooNE data, represented by the black points. The energy bins are 
numbered with the index j. The uncertainty is represented by the vertical error bars, which represent the sum of statistical and 
uncorrelated systematic uncertainties, (a) Expected number of !/e-like events N^fj calculated by the MiniBooNE collaboration. 



is given by the sum of the !/e-induced events {N^ 



= ,caK 



and the misidentified i^^-induced events [N^^'- ). (b) Best-fit 



JVC 

value of the number of i^e-like events Nl^°j obtained with the hypothesis of Ve disappearance. A^J|^j is given by the sum of 



N, 



,tho 

TabfH [MRv) 



UPl; 



U) 



M"^. 



and N.. 



J, the 



= f^N^ . The best-fit values of /^, sin 2-d and Am are those in the first column of 



[29l l30l |. PuJ-^u^^ is the survival probability of electron 
neutrinos in Eq. ([3]) averaged in the jth energy bin. The 
average in each bin is calculated using the ntuple-file of 
17,037 predicted muon-to-electron neutrino full transmu- 
tation events given in Ref. [29[, which contains informa- 
tion on reconstructed neutrino energy, true neutrino en- 
ergy, neutrino baseline and event weight for each event. 
The MiniBooNE measurement of a ratio 1.21 ± 0.24 of 
detected and predicted charged-current quasi-elastic v^ 
even ts Isif allows a value of f^ as large as about 15%. In 
Ref. [la] we used this estimate of the uncertainty of fi, in 
order to constrain its value in the least-squares analysis. 
Here we use directly the v^ data given in Ref. [231 ^^^ ^^^ 
construction of the MiniBooNE least-squares function 
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Xmb 



.-E 









J = l 



E K?. - ^3) (y^-:h>^ {K% - K7.) 



j,k=i 



(10) 



Here N°^^; are the numbers of measured i/p-like events 
in 11 reconstructed neutrino energy bins and N°^^- are 
the numbers of measured v^ charged-current quasi-elastic 



events in 8 reconstructed neutrino energy bins. The the- 
oretical number of i'^ events in the jth energy bin is given 
by 



n: 



the 



f.K: 



cal 



(11) 



where N^^^ , is the number of f^ events calculated by 
the MiniBooNE collaboration [29]. In order to take into 
account the correct statistical uncertainty corresponding 
to the rescaling of the number of i'^ events due to /^ in 
Eq. (jll[) . we used the covariance matrix Vi, given by 



(Kj,. = (Kr).fe + (/- - 1) ^c/^-fc 



(12) 



covariance matrix of v^ events 



where V,':^^ is the 8 x 
presented by the MiniBooNE collaboration in Ref. [29|. 
We did not use the complete 19 x J.9 covariance matrix 
of Ve and v^j. events given in Ref. |29| because the corre- 
lations involving Vg events have been obtained without 
taking into account the energy-dependent disappearance 
of electron neutrinos that we want to test. Assuming 
the correlations given in that 19 x 19 covariance matrix 



would suppress the energy dependence of Pi 



U) 



There- 



fore, for the uncertainties CTj^^j in Eq. (|10p we used only 
the diagonal elements of the i^e covariance matrix V^^^ 
given in Ref. [2^, corrected by the change of statistical 



uncertainty corresponding to the variation of expected 



events due to f^, and Pi 



U) 



in Eqs. dZ])-®: 



(K^ 



cal\ 



N, 



the 



n: 



cal 



(13) 



cal 



jU^ ,cal 



with iv^5 = iv;;;;"' + iv;: 

The result of the minimization of Xmbi/ ^^ shown in 
Fig. ^p, in which the solid histogram corresponds to the 
best-fit values of f„ , sin^ 2-3 and Arn^ in the first col- 
umn of Tab. m From Fig. [Hd, one can see that the fit is 
acceptable for all the Ve energy bins, including the first 
three bins which are out-of-fit in Fig. [iji. In Tab. U we 
give separately the contribution to Xmbi^ °^ ^^'^ ^^^^ three 
low-energy v,, bins and the sum of the contributions of 
the other v,, energy bins and all the f^ energy bins. In 
this way one can see that with f^ = 1 and P^^^^^ = 1 
(Null Hypothesis), although the global value x^ = 19.7 
is compatible with the number of degrees of freedom, 
NDF = 19, almost all the x^ is due to the anomalous 
contribution 14.3 of the first three low-energy v^ bins, 
whereas the other 16 i^e and i'^ energy bins are overfit- 
ted, with the excessively small x^ contribution of 5.4. 
This overfitting, which is probably due to an overesti- 
mate of the uncertainties, remains in the fit of the data 
with our hypothesis of /^ > 1 and i^e disappearance. On 
the other hand, our hypothesis clearly explains the low- 
energy anomaly reducing the x^ contribution of the first 
three low-energy v^. bins to the acceptable best-fit value 
of 2.0. 

Figure [21 shows the allowed regions in the sin^ 2'd-Am^ 



plane and the marginal Ax 



X 



An 



's for sin 2-d 



and Am^, from which one can infer the corresponding 
uncorrelated allowed intervals. One can see that the in- 
dication in favor of neutrino oscillations is not strong, 
being at the level of about 79% C.L. {\.2a). It is inter- 
esting to notice that the best-fit value of Am? is about 
2 eV , which is approximately the same best-fit value ob- 
tained in Ref. [14!| from the fit of the neutrino data of 
Gallium radioactive source experiments and the antineu- 
trino data of the Bugey and Chooz reactor experiments 
under the hypothesis of v^ and j/g disappearance. The 
results of the combined analysis of MiniBooNE neutrino 
data and the data of these other experiments is discussed 
in the following Sections. 



GALLIUM RADIOACTIVE SOURCE 
EXPERIMENTS 



III. 



The GALLEX [11|, 123, |3i and SAGE [10|, |35H37l col- 
laborations tested the respective Gallium solar neutrino 
detectors in so-called " Gallium radioactive source exper- 
iments" which consist in the detection of electron neutri- 
nos produced by intense artificial ^^Cr and '^^Ar radioac- 
tive sources placed inside the detectors. Taking into ac- 
count the uncertainty of the cross section of the detection 
process Ve + ^^Ga — >■ ^^Ge -I- e~ estimated in Ref. ^q, 
the ratio R of measured and predicted ^^Ge event rates 
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sin 2^ 



Ax 



FIG. 2. Allowed regions in the sin^ 2i9-Am^ plane and 
marginal Ax^'s for sin'^ 2i? and Am^ obtained from the fit 
of MiniBooNE neutrino data. The best-fit point is indicated 
by a cross. 
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and the average ratio is 



n 



Ga 



0.86 



-fO.05 
-0.05 ■ 



(14) 
(15) 
(16) 
(17) 



(18) 



Thus, the number of measured events is about 2.7cr 
smaller than the prediction. 

The theoretical prediction of the rate is based on the 
calculation of the detection cross section presented in 
Ref. [33] ■ It is possible that a part of the observed 
deficit is due to an overestimation of this cross section 
[lOl [371 l39j , because only the cross section of the transi- 
tion from the ground state of ^-'^Ga to the ground state of 
^^Ge is known with precision from the measured rate of 
electron capture decay of ^^Ge to ^^Ga. Electron neutri- 
nos produced by ^-'^Cr and '^^Ar radioactive sources can 
be absorbed also through transitions from the ground 
state of ^^Ga to two excited states of ^^Ge, with cross 
sections which are inferred using a nuclear model from 
p -\- ^^Ga -^ ^^Ge + n measurements [40|. This calcula- 
tion has large uncertainties [4l|, |42]. However, since the 
contribution of the transitions to the two excited states 



MBi/ Ga MBi/+Ga Re+^H (MB:/+Ga)+(Re+^H) 



Null Hyp. 


x' 


14.3 + 5.4 


9.4 




51.5 






NDF 


3 + 16 


4 




58 






GoF 


0.41 


0.051 




0.71 




Our Hyp. 


Amin 


2.0 + 7.6 


1.8 


2.2 + 9.2 


49.1 


4.1 + 63.4 




NDF 


16 
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20 


56 


78 




GoF 


0.89 


0.40 


0.93 


0.73 


0.80 




sin^ 2i9bf 


0.32 


0.27 


0.28 


0.042 


0.062 




Amgf 


1.84 


2.09 


1.92 


1.85 


1.85 




,/r 


1.26 




1.25 




1.17 


PG 


Ax4in 






0.098 


0.01 


6.97 




NDF 
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2 


2 




GoF 






0.95 


0.99 


0.03 



TABLE I. Values of x > number of degrees of freedom (NDF) and goodness-of-fit (GoF) for the fit of different combinations 
of MiniBooNE (MB;/), Gallium (Ga), and reactor (Re) data. The first three lines correspond to the case of /^ = 1 and no 
oscillations (Null Hyp.). The following six lines correspond to the case fi, > 1 and v^ disappearance (Our Hyp.). The last 
three lines give the parameter goodness-of-fit (PG) [32]. In the MSv column, the value of y^ in the Null. Hyp., the number 
of degrees of freedom in the Null. Hyp. (which is equal to the number of energy bins), and the value of Xmin in Our Hyp. 
are shown as the sum of the contributions of the first three low-energy Ve bins and the other v^ and i^^ energy bins. In the 
MBi/+Ga and (MBz/+Ga) + (Re+'^H) columns the value of Xmin in Our Hyp. is shown as the sum of the contribution of the 
first three MiniBooNE low-energy Ve bins and the other contributions. 



is only 5% [38|, even the complete absence of such tran- 



sitions would reduce Rp'^ to about 0.90lg;Q|, leaving an 
anomaly of about 1.8cr. 

Here we consider the electron neutrino disappearance 
explanation of the Gallium radioactive source experi- 
ments anomaly [iSl - KLSl l43l - t45| (another interesting ex- 
planation through quantum decoherence in neutrino os- 
cillations has been proposed in Ref. [2l[). 

In Ref. [14!] we have analyzed the data of the Gallium 
radioactive source experiments in terms of the effective 
survival probability in Eq. ([3]). Here we update that anal- 
ysis taking into account the revised value of Rq^^^^^ 
in Eq. (fT4|) published recently in Ref. [llj] and taking 
into account the asymmetric uncertainties of R^^^^^^^ 
Rq^^^^ and i?!^*^^ (which have been symmetrized for 
simplicity in the analysis presented in Ref. |14|). Fol- 
lowing the method described in Ref. [1J|, we obtained 
the best-fit values of sin^ 2i? and /S.rn? in the second col- 
umn of Tab.|I]and the allowed regions in the sin^ 2'd-lS.m? 
plane shown in Fig. [S] The indication in favor of neutrino 
oscillations is at the level of about 98% C.L. (2.3cr). 

From Tab. [2] and the comparison of Figs. [2] and [3] one 
can see that the fits of MiniBooNE and Gallium data 
lead to remarkably similar results: the best-fit values of 
the oscillation parameters are very close and the allowed 
regions in the sin 2'!?-A?7i^ plane are highly compatible. 
This is certainly an impressive success of our hypothesis 
of electron neutrino disappearance. 

The results of the combined fit of MiniBooNE and Gal- 
lium data are shown in the third column of Tab. |T] and 
in Fig. m The separate data sets are well fitted by the 
electron neutrino disappearance hypothesis: the y^ con- 
tribution of the first three MiniBooNE low-energy v^ bins 
is 2.2, that of the other 16 MiniBooNE v^ and v^ energy 
bins is 7.4, and that of the 4 Gallium data is 1.9. The 
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FIG. 3. Allowed regions in the sin^ 2i9-Am^ plane and 
marginal Ax^'s for sin^ 2i} and A.m? obtained from the com- 
bined fit of the results of the two GALLEX ^^ Cr radioactive 
source experiments and the SAGE ^^Cr and Ar radioactive 
source experiments. The best-fit point corresponding to Xmin 
is indicated by a cross. 



consistency of the combined fit is also supported by the 
excellent value of the parameter goodness-of-fit. Com- 
bining the two data sets improves the indication in favor 
of neutrino oscillations to the level of about 99.5% C.L. 
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FIG. 4. Allowed regions in the sin^ 2i9-Am^ plane and 
marginal Ax^'s for sin^ 2i? and Am^ obtained from the com- 
bined fit of the results of MiniBooNE neutrino data and the 
data of the gallium radioactive source experiments. The best- 
fit point corresponding to Xmin is indicated by a cross. 



FIG. 5. Ax^ as a function of mp. The horizontal lines corre- 
spond to the indicated value of confidence level. The dashed 
and dotted lines have been obtained, respectively, from the 
results in Eqs. ((2l| and ([22} of the Mainz and Troitsk Tri- 
tium /3-decay experiments. The solid line is the result of the 
combined fit. 



(2.8cr). 



IV. REACTOR AND TRITIUM EXPERIMENTS 

The indication of electron neutrino disappearance that 
we have found from the analysis of MiniBooNE and Gal- 
lium data must be confronted with the results of reactor 
electron antineutrino experiments. Assuming CPT in- 
variance, the survival probabilities of neutrinos and an- 
tineutrinos are equal (see Ref. [8|]). Thus, we can combine 
directly the results presented in the previous Section with 
the results of the analysis of the data of the Bugey and 
Chooz reactor experiments obtained in Ref. [Tjl . We are 
encouraged in this task by the coincidence of the best-fit 
value of Am^ at about 2 eV^. 

In addition to reactor neutrino experiments, also Tri- 
tium /3-decay experiments give information on the masses 
and mixing of neutrinos through the measurement of the 
electron energy spectrum in the process 



^H 



'He 



(19) 



The most accurate measurements of the effective electron 
neutrino mass (see Refs. @,[al) 



mp 




(20) 



have been performed in the Mainz [46j and Troitsk [47|: 
m^ = -0.6±2.2±2.1eV^ (Mainz), (21) 



'/3 



2.3±2.5±2.0eV'' 



(Troitsk) . 



(22) 



These measurements can be interpreted and combined in 
order to derive upper bounds for the effective mass m^ 
through a x^ analysis in the physical region mi > 0. In 
Fig. [5] we plotted the corresponding Ax^'s as a function 
of rufj . One can see that 

mp < 2.3 eV (Mainz, 95% C.L.) , (23) 

mp < 2.0 eV (Troitsk, 95% C.L.) , (24) 

in approximate agreement with the corresponding values 
in Refs. [4^|43]- The combined upper bound is 

mp < 1.8 eV (Mainz+Troitsk, 95% C.L.) . (25) 

In 3+1 four-neutrino schemes, taking into account that 
the mass splittings among mi, m,2, m,^ in Eqs. ([T|) and 
([2|) are negligible for a measurement of mp at the scale 
of 0.1 - 1 eV, the effective mass is given by 
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In 3+1 schemes of the type in Eq. ([S]), we have 



mp > I C/e4 1 V Am^ 



(26) 



(27) 



The effective mixing angle in short-baseline electron neu- 
trino disappearance experiments is related to \Ue4\ by 
Eq. ([6]). Inverting this relation and taking into account 
that the value of \Ue4\ must be small in order to fit the 
data of solar neutrino experiments with neutrino oscilla- 
tions, we have 



\Ue. 



\'^ 



\J\ - sin^ 2?? 



(28) 



In Fig.|6]we show the limits in the sin^ 2d~lS.m? plane ob- 



tained from Eqs. d^T]) and (US]) and the resuhs in Eqs. (PT|) 
and ((22|) of the Mainz and Troitsk Tritium /3-decay ex- 
periments. Notice that for small values of sin^ 2-3 the 
bounds are practically linear in the log-log plot in Fig. [51 
because in this case |C/e4p ~ sin^ 2i9/4 and the inequality 
in Eq. (gT]) leads to 



log Am^ < 2 log 2 



21ogm;^^ 



logsin2 2i9, (29) 



where m^ is the upper bound for rap. 

In the fourth column of Tab. H] and in Fig. [6] we re- 
port the results of the analysis of the data of the Bugey 
and Chooz reactor experiments presented in Ref. |14lpl. 
with the addition in the analysis of the results of the 
Mainz and Troitsk Tritium /3-decay experiments, which 
affect the high-Am^ region. As already commented in 
Ref. [14|, the reactor data are compatible with both the 
Null Hypothesis of absence of electron antineutrino dis- 
appearance and Our Hypothesis of electron antineutrino 
disappearance, with a hint in favor of electron antineu- 
trino oscillations due to a Am^ of about 2 eV. 
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FIG. 6. Allowed regions in the sin^ 2i9-Am^ plane and 
marginal Ax'^'s for sin^ 2i9 and Am^ obtained from the com- 
bined fit of the results of the Bugey and Chooz reactor ex- 
periments and the results of the Mainz and Troitsk Tritium 
/3-decay experiments. The three lines in the upper-right cor- 
ner are the exclusion curves obtained from the results of the 
Mainz and Troitsk Tritium /3-decay experiments alone. The 
best-fit point corresponding to Xmin is indicated by a cross. 



V. COMBINED ANALYSIS 

The results of the combined analysis of MiniBooNE, 
Gallium, reactor and Tritium data are presented in the 
last column of Tab. H] and in Fig. [71 One can see that 
the goodness-of-fit is high. The separate data sets are 
fitted fairly well by the electron neutrino disappearance 
hypothesis: the y^ contribution of the first three Mini- 
BooNE low-energy v^ bins is 4.1, that of the other 16 
MiniBooNE v^^ and i/^ energy bins is 7.5, that of the 4 
Gallium data is 6.3, that of the 56 reactor degrees of free- 
dom is 49.0 and that of the 2 Tritium degrees of freedom 
is 0.57. On the other hand, the 3% parameter goodness- 
of-fit of the combined analysis of neutrino MiniBooNE 
and Gallium data and antineutrino reactor and Tritium 
data is rather low. 



^ As erratum, let us notice that in tiic fourtli column of Table III in 
Ref. 11411 there is a small mistake in the evaluation of the param- 
eter croodness-of-fit. The correct values are Ay?,. = 0.52 and 
GoF = 0.47. We also notice that in the version of Ref. 11411 pub- 
lished in Phys. Rev. D the value of sin'^ 2i9hf for the Ga-|-Bu-|-Ch 
analysis (last column of Table III in Ref. [l4|) is different from 
the correct one, which is 0.054 (see the arXiv version of Ref. [ijl)- 



This low compatibility of the neutrino and antineutrino 
data sets is illustrated in Fig. |51 where we have plotted 
the marginal Ax^ 's for sin^ 2-3 obtained with the analysis 
of different data setfl One can see that 

1. The neutrino MiniBooNE and Gallium data agree 
to indicate a value of sin^ 2^ between about 0.11 
and 0.48 at 2a. 

2. The antineutrino reactor data indicate a value of 
sin^ 2^ smaller than about 0.10 at 2a. The Tritium 
data are practically irrelevant for the determination 
of sin^ 2d. 

3. The combined analysis is dominated by the reactor 
data and indicates a value of sin 2i? between about 
0.01 and 0.13 at 2a. 

The discrepancy between the neutrino and antineutrino 
determinations of sin^ 219 is about 2cr, in rough agreement 
with the above-mentioned 3% parameter goodness-of-fit 
of the combined analysis. In fact, the 2a disagreement 



^ We thank the anonymous referee of Phys. Rev. D for suggesting 
this interesting figure. 



- 68.27% C.L. {1a) 
• 95.45% C.L. (2a) 
•• 99.73% C.L. (3a) 
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FIG. 7. Allowed regions in the sin^ 2i9-Am^ plane and 
marginal Ax^'s for sin^ 2i9 and Am^ obtained from the com- 
bined fit of the results of MiniBooNE, Gallium, reactor and 
Tritium experiments. The best-fit point corresponding to 
Xmin is indicated by a cross. The three lines in the upper- 
right corner give the la, 2a and 3a limits in the sin^ 2{)-A.w? 
plane obtained from Eqs. (|20|l and (|28|) and the results in 
Eqs. ((21]) and ((22]) of the Mainz and Troitsk Tritium ^-decay 
experiments. 



between the neutrino and antineutrino data sets is en- 
tirely due to the different requirements on the value of 
sin^ 2z?, whereas they nicely agree on a best- fit value of 
Am^ at about 2eV^. 

Since the 3% parameter goodness-of-fit of the com- 
bined analysis shows a tension between the neutrino and 
antineutrino data (under our i^e-disappearance hypothe- 
sis) but is not sufficiently small to reject with confidence 
the compatibility of the neutrino and antineutrino data 
setQ in the following part of this Section and in Sec- 
tion I VII we consider the results and implications of the 
combined analysis. In Section IVIII we consider a possi- 
ble difference between the effective mixing angles in the 
neutrino and antineutrino sectors. 

Although the combined analysis of neutrino and an- 
tineutrino data favors smaller values of sin 2^ than those 
obtained from the analysis of MiniBooNE and Gallium 
data alone, the fit of the MiniBooNE and Gallium data 



FIG. 8. Marginal Ax^ = x^ 



Xn 



for sin 2i9 obtained from 



For example, the review on Statistics in the 2000 edition of 
the Review of Particle Physics [4g| says that if the goodness-of- 
fit "is larger than an agreed-upon value (0.001, 0.01, or 0.05 are 
common choices) , the data are consistent with the assumptions" . 



the analysis of different combinations of MiniBooNE, Gal- 
lium, reactor and Tritium data. The marginal Ax^ obtained 
from the analysis of reactor data alone and that obtained from 
the combined analysis of reactor and Tritium data are shown 
by the same line, since they practically coincide. 



remains better than in the case of no oscillations and 
/. = 1. 

Figure [H] shows the fit of MiniBooNE v^ data corre- 
sponding to the best-fit result of the combined analysis. 
One can see that the fit of the first three low-energy bins 
is not as good as that in Fig. [T]d, but it is nevertheless 
acceptable and much better than that in Fig. [1^. 

For the Gallium source experiments, the best-fit val- 
ues of the oscillation parameters give /jGALlex _ 



dGALLEX 
^Cr2 



0.97, i?g^"^ = 0.96 and R\^^^ = 0.96. 
Therefore, the experimental values of Rq^^^^ and 
i?c^°^ in Eqs. ^^ and ([16]) are fitted very well and 
the Gallium y^ contribution of 6.3 is almost equally due 
to the loose fits of R'^^^^^^ and i?|f ^^ in Eqs. US]) and 

Considering the combined fit of the results of Mini- 
BooNE, Gallium, reactor and Tritium data as a fair indi- 
cation in favor of a possible short-baseline electron neu- 
trino disappearance generated by the effective mixing pa- 



rameters Ar 



2eV and 0.01 < sin^ 2-d < 0.13, in the 



next Section we present the corresponding predictions for 
the effective neutrino masses in /3-decay and neutrinoless 
double-/3-decay experiments which could be measured in 
future experiments. 
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FIG. 9. Expected number of MiniBooNE Vi, events in the 
best-fit result of the combined analysis of MiniBooNE, Gal- 
lium, reactor and Tritium data (last column in Tab.|l]). The 
notation is the same as in Fig. [l] 



VI. PREDICTIONS FOR BETA-DECAY AND 

NEUTRINOLESS DOUBLE-BETA-DECAY 

EXPERIMENTS 

In this Section we present predictions for the effec- 
tive neutrino masses in /3-decay and neutrinoless double- 
/3-decay experiments obtained as a consequence of the 
combined fit of MiniBooNE, Gallium, reactor and Tri- 
tium data discussed in the previous Section. 



Figure [TO] shows the residual Ax^ 



X 



Cmin 



function of the contribution \UeA\\^rr? to the effective 
mass m^ in /3-decay experiments (see Eq. ([?7)) ). Since 
from the last column of Tab. Uwe have sin^ 2i?bf <C 1, we 
obtain 



\U, 



sin-^ 2^ 



bf 



e4|bf ■ 



= 0.016, 



and the best-fit value of |C/e4|vAm^ is 

(|C/e4|/ 

and 



Am2) =0.17eV, 

bf 



0.06 < |C/e4|VAm2 < 0.49 eV at 2cr . 



(30) 



(31) 



(32) 



This prediction is relevant for the KATRIN experiment 
|49| . which is under construction and scheduled to start 
in 2012. The expected sensitivity of about 0.2 eV at 
90% C.L. may be sufficient to observe a positive effect 
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FIG. 10. Ax" = x" 



Xmin 



as a function of the contribution 



I C/e4 1 \/ Ar7i2 to the effective /3-decay electron-neutrino mass 
m/3 in four-neutrino schemes obtained from the analysis of 
MiniBooNE and Gallium data (dashed line), from the anal- 
ysis of reactor and Tritium data (dotted line) , and from the 
combined analysis of the two sets of data (solid line). 



if |[/e4| vArr? is sufficiently large, as allowed by Ax^ in 
Fig.[TOl 

If massive neutrinos are Majorana particles, neutrino- 
less double-/? decay is possible, with a decay rate pro por - 
tional to the effective Majorana mass (see Refs. |2|. la. Isol- 

il) 



m2i3 



k 



Uek-mk 



(33) 



The results of the combined fit of MiniBooNE, Gallium, 
reactor and Tritium data discussed in Section |IV] allow 
us to estimate the contribution of the heaviest massive 
neutrino 1/4 to m2p, which is approximately given by 
|[/e4p vAtt?, taking into account the mass hierarchy in 
Eq. dSJ. 



Figure [TT] shows Ax^ 



An 



as a function of 



the contribution |C/e4pvAm2 in four-neutrino schemes 
to m2/3. The best-fit value is: 



C/e4|VAm2 =0.02eV, 



hi 



(34) 



and 



0.003 < |C/e4|VAm2 < 0.07 eV at2cr. (35) 



This range must be confronted with the expected con- 
tributions to ?7i2^ coming from the three light massive 
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neutrinos vi, V2, v^- Assuming a hierarchy of masses, 

nil ^ "T-2 ^ TTT-'i ^ "THi , (36) 

which is the most natural case compatible with the hier- 
archy in Eq. ([5]), we have 



^2^ 



t^e2VAm2oL 



ui^d^ml^^ + ul^^'^^ 



(37) 

where we have neglected the contribution of the lightest 
massive neutrino vi. From the 3cr upper limits of the 
three-neutrino mixing parameters given in Ref. [4|, we 
obtain 



2 X 10"-^ < It/, 



e2| 



'SOL 



<4x lO^-'cV, 



\U, 



e3| 



'ATM 



< 3 X 10"^ eV . 



(38) 



(39) 



Therefore, strong cancellations between the contributions 
of V2 and v^ are possible (albeit not likely [531), whereas 
the range in Eq. (|35|) disfavors strong cancellations be- 
tween the contributions of 1^2 and 1/3 and the contribu- 
tion of v^. In this case, m2^ ~ |C/e4pvAm^ leading to 
a possible observation of neutrinoless double-/? decay in 
future experiments which will be sensitive to values of 
7712,3 smaller that IQ-ieV (e.g. CUORE 54], EXO [H, 
SuperNEMO jll]; see the review in Ref. ^]). 

On the other hand, if neutrinoless double-/? decay ex- 
periments which are sensitive to values of m2p of the 
order of IQ-ieV (e.g. CUORICINO 1^, GERDA [El, 
Majorana [53|; see the review in Ref. [53) will see a pos- 
itive signal, maybe compatible with the signal asserted 
in Ref. \^Q\, the mass hierarchy in Eq. (p6| will become 
unlikely and the favorite 3-1-1 four-neutrino schemes will 
be those in which the three light neutrinos vi , 1^2 and v^ 
are almost degenerate at the mass scale of 7722/3. 



VII. MIXING ANGLE ASYMMETRY? 

The tension between neutrino and antineutrino data 
discussed in Section |V] could be due to a difference of the 
effective mixing angles in the neutrino and antineutrino 
sectors. Such a difference could be due to a violation 
of the fundamental CPT symmetry or to another un- 
known mechanism. Phenomenological analyses of differ- 
ent masses and mixings for neutrinos and antineutrinos 
have been presented in several publications J6ll - I-74| . 

In this section we consider the possibility that neutri- 
nos and antineutrinos have different effective masses and 
mixings in short-baseline v^. and i>e disappearance exper- 
iments. We fit the neutrino and antineutrino data with 
the survival probabilities 



P^^^^^ iL,E) = l- sin^ 2i?, sn 



Pp°Xp, iL,E)^l~ sin^ 27?p sh 



AmlL 

AE 
AmlL 

AE 



(40) 
(41) 
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FIG. 11. Ax^ = X'^ ~ Xmin ^s a function of the contribu- 
tion |(7e4p\/Am2 to the effective neutrinoless double-/? decay 
Majorana mass 7712/3 in four-neutrino schemes obtained from 
the analysis of MiniBooNE and Gallium data (dashed line), 
from the analysis of reactor and Tritium data (dotted line), 
and from the combined analysis of the two sets of data (solid 
line) . 



The results for the two fits are those presented in Fig. |4] 
and the third column in Tab. U for neutrinos (MBj^+Ga) 
and Fig. [6] and the fourth column in Tab. U for antineu- 
trinos (Re-t--^!!). 

Since the fit of the data does not require a difference 
of Ami and Am?, we consider only the mixing angle 
asymmetry 



A, 



'■2§ 



sin^ 2i% - sin^ 27?p . 



(42) 



Figure [12] shows the marginal Ax^ as a function of 
v4si„2 2ij. The best- fit value of ^sin^ 21? is 



Ain2 2i5 — 0-23: 



and the 2cr allowed range of A^^^2 2^ is 
0.06 < ^i„2 2^ < 0.45 



(43) 



(44) 



but there is no limit on the asymmetry at 3a. The sta- 
tistical significance of ^31^2 2^ > is 99.14% C.L. (2.6a). 
It is interesting to note that a difference between neu- 
trino and antineutrino mixings can be tested in /3-decay 
experiments by searching for different effective neutrino 
masses in /3^ and /?+ decays. The prediction for the con- 
tribution ( I Ue4 1 vAm? j to the effective electron an- 
tineutrino mass in /3~ decays from the analysis of an- 
tineutrino reactor and Tritium data can be obtained from 
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FIG. 12. Marginal Ax^ = X^ — Xmin as a function of the 
mixing angle asymmetry A^'^^2 2^. 



the dotted line in Fig. [101 the best fit is 



{yjeAW^rnA =0.14eV, 



(45) 



and 



{\UeiWl^mA <1.29eV at 2cr . (46) 



to the effective 



For the contribution ( | [/e4 1 vAn 

electron neutrino mass in /3+ decays we must consider 
the dashed line in Fig. I10[ which has been obtained from 
the analysis of MiniBooNE and Gallium neutrino data: 
the best fit is 



(I 



bf 



[/e4|%/A^i2) =0.38eV, 



and 



(\Uei\^/'^^) ^ 



>0.21eV at2cr. 



(47) 



(48) 



Unfortunately the existing and foreseen experiments 
are /3~-decay experiments for which the contribution 

|C/e4|vAm^) to the effective electron antineutrino 



(|C/e4|yA^)^ 



mass is expected to be small. The future /3~-decay ex- 
periment will use either Tritium (KATRIN [49|) or ^^^Re 
(MARE [73 ) ■ If the mixing difference between the neu- 
trino and antineutrino sectors will be confirmed with 
highest confidence by future neutrino oscillation data it 



will be interesting to study the possibility of making /3+- 
decay experiments for the search of the effective electron 
antineutrino mass, for which the dashed line in Fig. [10] 
and Eq. (pS)) give a reachable lower limit. 

Wc do not consider here neutrinoless double-/3 decay 
in the case of a neutrino-antineutrino mixing difference, 
since the Majorana nature of neutrinos requires a treat- 
ment whi ch g oes well beyond the purposes of this paper 
(see Ref. [76f). 



VIII. CONCLUSIONS 

In this paper we have discussed a neutrino oscillation 
interpretation of the MiniBooNE low-energy anomaly 
and the Gallium radioactive source experiments anomaly 
in the framework of 3-1-1 four- neutrino mixing schemes. 
We have shown that the combined fit of MiniBooNE and 
Gallium data indicate a possible short-baseline electron 
neutrino disappearance generated by effective oscillation 
parameters Am^ > 0.1 eV^ and 0.11 < sin^ 2-d < 0.48 at 
2cr, with best fit at Am^ ~ 2eV^ and sin^ 2i) ~ 0.3 (see 
Fig. HI). 

We have also considered the data of the Bugey and 
Ghooz reactor neutrino oscillation experiments and the 
results of the Mainz and Troitsk Tritium /3-decay experi- 
ments, which imply an upper bound on the effective elec- 
tron neutrino mass of about 2 eV (see Fig. [S] and the com- 
bined upper bound in Eq. (US])). As already discussed in 
Ref. [ij , the Bugey data give a faint indication of a pos- 
sible short-baseline electron neutrino disappearance gen- 
erated by effective oscillation parameters Am^ ~ 2 eV^ 
and sin^ 2i5 ~ 0.04, which is compatible with Chooz and 
Tritium data (see Fig. [6]). 

In Section [V] we have discussed the tension between 
the neutrino MiniBooNE and Gallium data and the an- 
tineutrino reactor and Tritium data. Considering such 
tension as a statistical fluctuation, we have presented 
the results of the combined analysis of MiniBooNE, Gal- 
lium, reactor and Tritium data: Am^ ~ 2 eV and 
0.01 < sin^ 21} < 0.13 at 2cr, with best fit at Am^ ~ 2 eV^ 
and sin^ 2z9 ~ 0.06 (see Fig. [7]). 

In Section Jv!] we have presented predictions for the 
effective neutrino masses in /3-decay and neutrinoless 
double-/3-decay experiments obtained as a consequence 
of the combined analysis of MiniBooNE, Gallium, reac- 
tor and Tritium data, assuming the hierarchy of masses 
in Eq. ^ . The predicted interval for the contribution of 
77i4 to the effective neutrino mass in /3-decay is between 
about 0.06 and 0.49 eV at 2a. The upper part of this in- 
terval may be reached by the KATRIN experiment |49{ . 
For neutrinoless double- /3-decay we obtained a prediction 
for the contribution of 7714 to the effective neutrino mass 
between about 0.003 and 0.07 eV at 2a which may be 
reached in future experiments (see Ref. [52|). 

We also considered, in Section IVIIl the possibility of 
reconciling the tension between the neutrino MiniBooNE 
and Gallium data and the antineutrino reactor and Tri- 
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tium data discussed in Section|V]'with different mixings in 
the neutrino and antineutrino sectors. We found a 2.6a 
indication of a mixing angle asymmetry (99.14% C.L.). 
We pointed out the possibiHty of checking the mixing 
difference between the neutrino and antineutrino sectors 
by measuring different effective electron antineutrino and 
neutrino masses in /3~ and /3+ decay experiments. 

The indication in favor of short-baseline disappearance 
of electron neutrinos imply the possible existence of a 
light sterile neutrino which could have important conse- 
quences in physics [l8l42q . |22| . I77H83J , astrophysics [8J- 
[93 | and cosmology |93l - [98j . 

As far as the effective number of neutrino species in 
cosmology, A^cff, is concerned, the analysis of 7-years 
WMAP data has provided the following result: iVoff' = 
4.34to;88 (68% C.L.) 0. In 2011 the Planck experi- 
ment will measure A'cft with a factor of 4 improvem ent 
in accuracy with respect to present data |lOOl Il0l| . In 
other words, the possibility of existence of a fourth light 
sterile neutrino could be pursued with 5a significance. 

Finally, we would like to encourage all experiments 
which can investigate the hypothesis of short-baseline 
electron neutrino disappearance. 

Starting from 2010, at the same L/E of MiniBoone, the 
magnetic off- a xis near detector at 280 m of the T2K ex- 
periment [l0.2J will count fe events with expected higher 
statistics and similar v^ background contamination. A 
test of short-baseline oscillations may be done, although 
the accuracy suffers from the scarce knowledge of the 
neutrino flu x and of the neutrino cross section at 1 GeV 
energies [l5lll03| . 

A better measurement will be po ssible with the new 
CERN-PS neutrino beam flol[lOa|. thanks to the pres- 
ence of 2 detectors at 140 m (NEAR) and 885 m (FAR). 



At Af 



2 eV , the oscillation length is about 1 km 



for 1 GeV neutrino energies. Therefore, one can reduce 
the systematic error of the Monte Carlo predictions by 
normalizing the high energy part of the v^ spectrum at 
the NEAR location. In addition, a better j/^ background 
rejection will be possible using the liquid Argon technol- 
ogy. The interesting possibility of a i^e tagging in the 
CERN-PS beam was also studied in this context [106]. 

New measurements with a radioactive source could be 
made in the SAGE experiment [Ifll, with the Borexino 
detector and with the future LENS detector [73]. At 
Am^ a: 2eV^, the oscillation length is about 1 m for 1 
MeV neutrino energies. Therefore Borexino could mea- 
sure the oscillation pattern over a distance of 4 m (the 
Borexino radius) using the well known v^-e scattering 
process and wit h a v ertex resolution that at the moment 
is about 15 cm jl07l |. 

At the Gran Sasso laboratories a very interesting mea- 
surement could be realized by using the ICARUS 600 
ton detector and new low-cost and high-power proton 



cyclotrons under development for commercial uses [lOSl ] . 
These provide electron neutrino beams with energy up 
to 52 MeV from muon decay- at-rest. A low-energy I'e 

disappearance experiment (as well as f^ disappearance 
and ly^ — >■ u^ measurements) can be performed with such 
devices due to the full efhciency of the ICARUS detec- 
tor at 20 MeV energies. The expected event rate is about 
400 charged-current electron neutrino events per year per 
ton wit h th e ICARUS detector located at 50 m from the 
source |l09l | . The number of events is calculated assum- 
ing 10^^ lye^s per year and a fully efhcient detector. Since 
at Am^ ~ 2 eV^ the oscillation length is about 20 m for 
a 20 MeV neutrino energy, it is possible to measure the 
full oscillation pattern along the beam direction inside 
the ICARUS volume. 

The disappearance of electron neutrinos can be inves- 
tigate d wit h high accuracy in future nea r-detector beta- 
beam jllOJ I and neutrino factory [73, llll| experiments in 
which the neutrino fluxes will be known with high preci- 
sion. 

Furthermore, the MiniBooNE low-energy anomaly 
may be clarified by the ArgoNeuT, MicroBooNE [11J| 
and BooNE experiments [113J | , and the magnetic off-axis 
near detector at 280 m of the T2K experiment has the 
unique opportunity to mea sure the charge of the events 
of the low-energy anomaly J114l | . 



NOTE ADDED 

After the completion of this work, two important ex- 
perimental results have been presented at the Neutrino 
2010 conference: 

1. The MiniBooNE collaboration presented updated 
results on the search for short-baseline Df^ — >■ Dp 
oscilla tions which are compatible with the LSND 
signal J115l Ill6| . The inclusion of these data in our 
framework will require a separate analysis in which 
the assumption of negligible |C/^4p is relaxed [117J . 

2. The MINOS collaboration presented an indication 
of a possible difference between the effective mix- 
ings of neutrinos and anti neutr inos in long-baseline 



f^ and i'^ disappearance |118l |. This indication is 
analogous to that discussed in Section FVlII 
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